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Abstract. We study the 3D Edwards-Anderson model with binary interactions 
by Monte Carlo simulations. Direct evidence of finite-size scaling is provided, and 
the universal finite-size scaling functions are determined. Using an iterative ex- 
trapolation procedure, Monte Carlo data are extrapolated to infinite volume up 
to correlation length ^ « 140. The infinite volume data are consistent with both 
a continuous phase transition at finite temperature and an essential singularity at 
finite temperature. An essential singularity at zero temperature is excluded. 

1 Introduction 

Understanding the thermodynamics of spin glasses represents a challenging 
unsolved problem for statistical and computational physics |Q] . The very exi- 
stence of a phase transition in the 3D Ising spin glass is still unclear. Previous 
Monte Carlo (MC) simulations |^,H,^J^,^ give a certain evidence of a Tc 7^ 
continuous transition with power-law divergence of ^ at T ^ . However, 
they cannot exclude completely neither an exponential divergence at T — s- 0^ 
, nor a line of critical points at T < Tc ^ §,§,01 , which implies an 
exponential divergence at T — > T^. Furthermore, the important topic of 
finite-size scaling (FSS) has not been much investigated for spin glasses. 

Here, we study the 31? Ising spin glass by MC simulations of moderate- 
size systems in the paramagnetic phase (less hampered by long equilibration 
times), using the powerful FSS method of In Sect. 2 details on the sim- 
ulations are given. In Sect. 3 a direct evidence of FSS, independent of the 
nature of the divergence of ^, is provided, and the universal FSS functions 
are determined for the first time [^. The MC data are extrapolated to infi- 
nite volume up to ^ one order of magnitude larger than before. In Sect. 4 the 
critical behavior of the extrapolated data is analyzed. 

2 Model and simulation 

We consider the 3D Edwards- Anderson model, whose Hamiltonian is 

H = — ''^a^JxyCTy (1) 

where ax are Ising spins on a simple cubic lattice of linear size L with periodic 
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boundaries, and J^y are independent random interactions taking the values 
±1 with equal probability. The sum runs over pairs of nearest neighbor sites. 
We simulate with the heath-bath algorithm many samples of the model 
with different J^y- From two independent replicas (cr, t) with the same 
Jxyi we measure the overlap qx = OxTx and q = L^^ qx- We compute the 
second-moment correlation length using the following definition: 



V 4 sin2(p/2) 

where S{k) is the Fourier transform of the overlap correlation function 

^(fc) = ^e*'^(q,g,+.) , (3) 

r 

(the arguments T, L are omitted) and p = (27r/L,0,0) is the smallest non- 
zero wave vector. We also measure the spin-glass susceptibility xsg{T, L) = 
L^{q'^) — S{0). The runs are done on a Cray T3E with a fast code (see 
[|j for details) that exploits the high parallelism of spin glass simulations, 
combining multi-processor parallelism with an efficient multi-spin coding 
technique. Average speed on a single processor (DEC Alpha EV5, 600 MHz) 
is 4.5 X 10'' spin updates per second, and we typically used 32 to 128 PEs. 
We simulated 104 pairs (T, L), (T, 2L) with L between 4 and 48. 



3 Finite size scaling analysis 

According to the FSS hypothesis [|lO), if 0{T,L) is some long-distance ob- 
servable (for example, £,{T,L) or xsg{T, L)), then 



0{T, oo) 



/o(e(r,^)/i) (4) 



where fo is a universal function and corrections to FSS are neglected. From 
(^ we obtain a useful relation involving only finite-volume observables 

°(^-^« . FJiiTMIL) (5) 



OiT,L) 

where Fq is another universal function. As shown in Fig|^a,b and Fig.^, our 
data for the observables XsGi <1a = (9^) and F = V^^S{p) verify the 
ansatz (||). This provides a direct test of the FSS hypothesis. Small system- 
atic deviations, due to corrections to FSS, appear in Fig.^ for L = 5 and 
^{T,L)/L = X > 0.4. We emphasize that FSS is not assumed a priori and 
that no adjustable parameters are contained in (H). Furthermore, no partic- 
ular dependence of the observables on the temperature is assumed. 

We fit the data in Fig.|l|a,b to two suitable functions F^^^ , and then 
extrapolate the pair {S,,xsg) from L ^ 2L ^ 2^L . . . ^ oo vising (^ 
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Fig. 1. (a,b) FSS plot with the form in (||) for (a) O = XSG and (b) O = Error 
bars (estimated with ajackknife procedure) are one standard deviation. (c,d) FSS 
plot with the form in for (c) O = xsG and (d) O = ^. See for how to estimate 
error bars of the extrapolated data. 



iteratively. For any temperature, we verify that extrapolations from different 
L agree within the error bars. In this procedure, we assumed implicitly that 
with a given function Fq will continue to hold as L oo. This assumption 
could fail if there is a crossover at large L. However, at high T extrapolations 
from small L are consistent with data from large L, which have almost no 
finite-size effects, and thus a crossover is unlikely. We have a good control on 
the extrapolated data up to ^ ~ 140; at lower temperatures the statistical 
errors become quite large, and the data are more sensitive to FSS corrections. 

In Fig.|l|c,d we show that with our extrapolated data (|^) is satisfied re- 
markably well, providing a further test of the method. If O ^ ^10 /f ^ _> 
then fo{x) in (^) must satisfy fo{x) ~ x~'^°^^ as a; — > 00. As shown in the 
insets of Fig.|l| (c,d), our curves indeed have a power-law asymptotic decay. 
We emphasize that all the scaling functions in FigsJI^ and are universal. 
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Fig. 2. FSS plot with tlie form in (|) for (a) O = qi = (g"*) and (h) O = F 
V-''S{p), \p\ = 2n/L. 



4 Nature of the phase transition 

We now compare our extrapolated data with the three different scenarios 
compatible with previous simulations. 

(i) Power-law .singularity at Tc =/= 0. We fit our data to 

e(r) = c^T - T,)-^ [1 + (T - nf] (6) 
XsaiO^be-'^ [l + dC^] (7) 

with fixed 9 and A. Without corrections to scaHng (a^ ~ d = 0), the fit 
parameters show small systematic variations when we vary the fit interval. 
The fits stabilize with 1 < 6* < 2 and 1 < Zi < 1.5, the preferred values being 
9 = 1.4 and A = 1.3 {goodness of fit parameter Q > 0.6 and Q > 0.98 resp.). 
Our estimates for the critical parameters are Tc — 1.156±0.G15, i' = 1.8±0.2, 
and r] = —0.26 ±0.04. Corrections to scaling are important for f < 10 (Fig.||). 

(ii) Essential singularity at Tc ^ 0. Our data fit very well also to 

e(r) = exp(g4/(T - T^r) . (8) 

For ^ > 3.8 the best fit (shown in Fig.|b) gives a = 0.5 ±0.3, = 1.08 ±0.04 
{Q = 0.69). Deviations from this fit for f < 3 are consistent with corrections 
to scaling of « 10%. Since for an exponential singularity we expect multi- 
plicative logarithmic corrections, we tried also the fit 

XSGiO = bi e-'^' {logO^ (9) 

obtaining r]i = -0.36 ± 0.03 and r = -0.36 ± 0.06 (Q > 0.9) (see Fig.|a). 

(iii) Essential singularity at T = 0. Fitting our data to 



e(r) = /4 exp(g4/T-) 



(10) 
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we find a > 9. Such a high value is highly implausible on the basis of renor- 
malization group arguments, from which we expect cr « 2 Q. We therefore 
believe that an essential singularity at T = is excluded. 
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Fig. 3. Critical behavior of the infinite volume data, (a) Best fit to ([^) for ^ > 1.8 
(line A), leading term from the same fit (A leading), best fit to (^) for ^ > 2.2 (B). 
(b) Best fit to (^) for ^ > 1.9 (A), leadin g te rm from the same fit (A leading), best 
fit to d) for ^ > 3.8 (B), and best fit to (|l^ for ^ > 14 (C). 
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